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Abstract 

The order a 2 s contribution to the inclusive semileptonic decay width of a b 
quark T(b — > X u eV e ) is calculated analytically for zero mass u quarks 
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Measurements of the semileptonic B meson decay rates at B factories [1[] provide 
precise means to extract the Cabibbo-Kobayashi-Maskawa matrix elements \V u b\ and 
\V c b\ from experiment. The accurate determination of these parameters is a requisite 
for many tests of CP violation within the three-generation Standard Model. 

The Heavy Quark Expansion provides a framework for the systematic calculation of 
contributions to the inclusive B meson decay rate [0-0. To leading order in the heavy 
quark mass the inclusive B meson decay rate is equal to the decay rate of a on-shell 
b quark treated within renormalization group improved perturbative QCD. The non- 
perturbative corrections are suppressed by at least two powers of the heavy quark mass 
and can be expressed as matrix elements of higher dimension operators in the heavy 
quark effective field theory. 

Perturbative QCD contributions to b quark decay play an important role in accurate 
predictions for the inclusive B meson decay rate. The first order QCD corrections to 
the inclusive semileptonic b — > u decay width Y(b — > X u eV e ) were derived from the 
calculation of the 1-loop QED corrections to the muon decay rate |7|]. For Y(b — > X c eV e ) 
c quark mass corrections are important and the 1-loop QCD corrections were obtained 
in a closed analytical form in Ref. ||. 

For dimensional reasons Y(b — > X u eV e ) is proportional to the b quark mass to the 
fifth power. The adopted definition for the b quark mass therefore strongly affects the 
coefficients of the higher order QCD corrections. If the width is expressed in terms 
of the MS renormalized b quark mass one finds that the 1-loop QCD corrections to 
Y(b — > X u eV e ) are significant, one-third of the size of the tree level contribution. 

The 2-loop QCD contribution to Y(b — > X u eV e ) and Y(b — > X c eV e ) proportional to 
the number of light quark flavours was calculated in Refs. P,[l0|. This result was used 
to obtain an estimate for the full 2-loop corrections for both semileptonic b —>■ c and 
b — > u decay ||10|| . Recently an accurate estimate of the full 2-loop correction to the 



inclusive semileptonic b — > c decay rate was obtained JTT[] based on calculations at 3 
distinct values for the invariant mass of the lepton pair [|TT| -|T3"|| . 

In this letter the 2-loop QCD corrections to the inclusive semileptonic b — > u decay 
rate are evaluated analytically. The starting point of the calculation is the Lagrangian 



C — C w + £qcd (1) 
Here Cw is the effective weak interaction term 

Q 

£w = --^V ub [ip u -f X (l - 75)^5] • [V^7a(1 - 7s)^ e ] (2) 

where ip u , ipb, 4>e and ip Ue are the wave functions for the u quark, b quark, electron 
and its neutrino respectively. £qcd is the standard QCD Lagrangian for the strong 
interactions. The calculation is performed in effective 5 flavour QCD since the t quark 
decouples for the process under consideration. The QCD correction to the semileptonic 
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Figure 1: Examples of non-abelian diagrams whose cuts give contributions to b 
b — ► ueV e g, b — > ueV e gg. 



ueiJ P 



decay rate are finite since they are separately finite for a vector quark current and the 
flavour non-singlet axial- vector quark current. 



Throughout this article we use dimensional regularization [|14| ,|15| and the standard 
modification of the minimal subtraction scheme [IB|, the MS-scheme [[T7[ for coupling 
constant renormalization. Calculations with massive external fermions are most easily 
performed in the on-shell renormalization scheme for fermion masses. Nevertheless, 
the final results in this letter will be expressed in terms of MS renormalized quark 
masses. For the treatment of the 75 matrix in dimensional regularization the technique 
described in Ref. |TB| is used which is based on the original definition of 75 in Ref. ||15|| . 
It is interesting to note that in the limit of massless u quarks the contribution to the 
decay rate coming from the axial-vector part of the weak interaction Eq. (0) can be 
reduced to the vector part by using the effective anticommutation property of the 75 
matrix. 

For the order a 2 s corrections to the inclusive width Y(b — > X u eV e ) one has to deter- 



mine the sum of the decay widths b 



ueu e , b 



ueu e g, b — > uev e gg and b — > ueu e qq, 



where g denotes a gluon and q denotes a quark, with up to two gluons or quarks in vir- 



tual corrections. According to the Kinoshita-Lee-Nauenberg [TjJ theorem, the inclusive 
decay rate is free from singularities as the u quark mass goes to zero. The calculation 
will be done in the approximation of massless u,d,s and c quarks. The c quark enters 
only in the partial rate b — > ueV e cc or as cc pairs in virtual corrections, and the effect 
of a non-zero c quark mass on the total decay rate is therefore limited. The situation is 
quite different for T(b — ► X c eV e ) where the c quark enters in all diagrams and c quark 
mass effects are very important. 

Using the optical theorem one may express the inclusive decay rate as the imaginary 
part of 4-loop on-shell propagator type diagrams. Some of these diagrams are shown in 
Fig 1. The imaginary part of these 4-loop quark propagator diagrams can be evaluated 
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analytically using the methods employed in Ref. [20J where the 2-loop QED corrections 
to the muon decay rate were evaluated. As was noted before, the muon decay rate and 
r(6 — > X u eV e ) are closely related; the result for the 2-loop QED corrections to the muon 
decay width corresponds formally to keeping only the abelian part in the 2-loop QCD 
corrections to F(b — ► X u eV e ). The approach of Ref. [20] relies heavily on the use of the 
method of integration-by-parts |21] within dimensional regularization to express the 
imaginary part of the 4-loop diagrams in terms of a small set of 26 primitive integrals. 
Although there are diagrams for F(b — > X u eV e ) that involve integration topologies that 
were absent in Ref. [^0J it turns out that after the use of integration-by-parts identities 
for these topologies no new primitive integrals are needed. 

All diagrams that contribute to T(b — > X u eV e ) are calculated in a general covariant 
gauge for the gluon fields. The explicit cancellation of the gauge dependence in the 
sum of the diagrams gives an important check of the result. The result that is obtained 
in this way reads 



T(b -> X u ev e 



G F \V ub \ 2 Ml 

1927T 3 



0(or 



Ti 



Ti 



(3) 



C F (f - 3C 2 ) 
-2.41307 



CaCf 



154927 
10368 

11047 



53 1 ]n(9\ 4- 9 V 38 V 4- 10 V 



275 j 
96 111 



M? 



^C 2 ln 



, r , r2 + 53C 2 ln(2) - ipC2 - ^Cs + f Ci] 



+C F T F n f [-ff + i C 2 + |C3 - | In (£) + C2 In ( 



+C F T F _ | Ca _ 24C 3 ) 



-21.29553 -4.625050 In 



where 5 (/V) is the coupling constant in effective 5 flavour QCD, /x is the renormal- 
ization scale, rif = 5 is the total number of quark flavours and Mb is the b quark pole 
mass in the two loop order. C F = 4/3 and Ca = 3 are the Casimir operators of the 
fundamental and adjoint representation of the colour group SU(3), T F = 1/2 is the 
trace normalization of the fundamental representation. denotes the Riemann zeta 
function, C2 = vr 2 /6, Cs = 1.2020569 • • • , C 4 = tt 4 /90. 

We note that the coefficient of C F T F nf in Eq.(U) agrees with the value 3.22 obtained 
in Ref. and the terms proportional to Cp and C F Tp agree with the result for the 
QED contributions to the muon decay width of Ref. [20|. The order a 2 s coefficient in 
Eq. (f|) is sizable but this coefficient is somewhat smaller than the estimate —28.7 
0- 



Mr 



(4) 



obtained in Ref. 
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It has been shown [0,|23| that the perturbative coefficients for T(b — > X u eV e ) receive 
contributions that grow rapidly at higher orders corresponding to a singularity in the 
Borel plane at n/ (2/3 ) (with /3 given below Eq. ([15])) when the decay width is expressed 
in terms of the b quark pole mass. This particular growth of the perturbative coefficients 
is related to the infrared sensitivity of the pole mass definition and is absent when the 
inclusive decay width is expressed in terms of the MS renormalized b quark mass mf,(/i). 
We will therefore express the decay width in terms of m^/i) using the known relation 



between the pole quark mass and the MS renormalized quark mass \2A, 25 



1 + ci — + c 2 | —) +0(af 



7T 



7T 



(5) 



ci = Cv(l + fm( 
CaCf 



c 2 



^ + §C 2 hi(2)-±C 2 -f<3 + Wm'^ 



384 



, , i I! l n 2 f Jft 

M b J + 32 



+cl 



^-3C 2 ln(2) + fC2 + IC3 + lln hg 



J 2S 



9 lr, 2 JL 



+ #ln 



+C F T F n f -gg - 5C2 - 2i In 



+ C F T F (-l + IC2) 



(6) 



Applying this relation to eliminate M& from Eq. (|3|) and from the logarithms appearing 
in Eqs. (|,D 



In 



JL. 

Ml 



In 



2C f 



ai 5) (//) 



3 

1 + - In . 

4 Vm 2 (/i)yj 7T 



(7) 



one obtains the result for the decay width in terms of the MS renormalized b quark 
mass 



, kl ^itl +h2 (^M\ + 0(a »> 



7T 



7T 



(8) 



5 



f -3C 2 + ^ln(^ 



m?(A«) 



4.25360 + 51n 



mf(/0 



C A C 



19057 



648 



-19C 2 ln(2) + |C 2 -lC3 + tC4 + fln(4 



2 

50 



.11 i,. / k 



+ lln 2 



™£(/0 



-fC 2 ln 
+C F T F rif 



281113 _i_ oo/- l Tl/ 'oN 9455a 22a , 67 a , 405 i„ f K 

10368" + ^ in(2) - w ^2 - ^Cs + + 7e" ln I Hn^ 



22", ,,,2 / ft 



+ ifln 



™|00 



1037 _ 13 A , 8 A _ 15 i 

36^2 T 3 <,3 4 V m g((u) y 1 V m b(^) 



144 36 



ft 
7TT 



+ C 2 ln 



fin 2 



mf(/0 



~\-CpTp (f§ + 3C 2 - 24C 3 ) 

26.78476 + 36.99016 In (^y) + 17.291671n 2 ( 



m£(/x) 



(9) 



where rrib{jj) = mf\fj) is the MS renormalized b quark mass in effective 5 flavour 
QCD. The decay width may also be expressed in terms of parameters of 4 flavour 
QCD by using the known decoupling relations for the coupling constant and the quark 
masses [p6| -[28[ 



7T 7T I 7T ) 3 UftW 



„( 4 ), 



+ o(o2) 



1 + 



«i 4) (/x) 



71 



1 



5 / /i 2 

H In ^ . 

24 V^(/^) 



89 
288 



T F C F ( --In 



mg(//) 



(10) 



to obtain 



«i 4) (/i) , / «i 4) (/i) 

1 + / 1 r J2 



71 



+ 0(a: 



71 



(12) 



/i 

h 



4.25360 + 5 In 



25.75467 + 38.39353 In 



m?(/i) 



+ 17.70833 In 2 



mf(A*) 



(13) 



Here m&(//) = rn£\n) and only the numerical expressions are given. 
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Each of the obtained expressions for the decay width Eqs. (I),© an d flT2| ) satisfies 
formal renormalization group invariance that is required of a physical quantity which 
in the a 2 approximation reads 



dln(/i 2 



-rV, a a (/x), m b (fj)) = 0(a 3 (/i)) (14) 



The fi dependence in the renormalized mass and coupling constants is given by the 
renormalization group equations 



da 



dlia.fi 
dham, 
din// 



= /3(a) = -f3 a 2 - ^a 3 - (3 2 a 4 - (3 3 a 5 + 0(a 6 ) (15) 
■= 7m(a) = -7o«-7ia 2 -72a 3 -73a 4 + 0(a 5 ) (16) 



where a = a s /n, [3 = ±±C A - \T F n f , 70 = §C F , 71 = C F (^C F + ^C A - ^T F n f ) and 
rif is the number of quark flavours in effective QCD. The higher order coefficients up 
to /3 3 and 73 are given in p9| . 

As is indicated in Eq. flT4]) the two loop QCD expressions for Y{b — > X u eV e ) have a 
remaining renormalization scale dependence that is of order a 3 . The renormalization 
scale dependence of the expression for the decay width Eq. (d) in the different orders of 
perturbative QCD is illustrated in Fig 2. In Fig. 2 the leading order curve (LO) corre- 
sponds to the tree level expression for the decay width and to solving the renormaliza- 
tion group equation for a s (/i) and m^fi) in the leading order, the next-to- leading order 
curve (NLO) corresponds to keeping one higher order in a s both in the expression for the 
decay width and the renormalization group equations, etc. The renormalization group 
equation for a s (fi) is solved numerically starting from af\mz) = 0.118. Similarly the 
scale evolution of mf\fi) is done numerically starting from m^m^) = 4.3 GeV. For 
comparison we have included a NNNLO curve that corresponds to taking the presently 
unknown constant of order a 3 in Eq. (R|) to be 5 3 = 125, assuming here an approximate 



power-like growth of the lower order coefficients, see e.g. Ref. |3(J. The logarithmic 
terms a 3 In* (fi 2 / 'ml(fi)) , i — 1, 2,3 can be expressed in terms of the known lower order 
coefficients by requiring formal /^-independence of the decay rate in the a 3 order. 

One may see from Fig. 2 that in the NNLO expression Eq. (|5p a moderate renor- 
malization scale dependence remains near /1 = m^mj). Variation of fi in the relatively 
wide interval between fi = 1/2 m^mi,) and fi = 2mb(rrib) changes the decay width in 
the NNL order by about 15%. If variation with respect to fi in this interval is used to 
estimate the contributions to T(b — >• X u eV e ) beyond the a 2 order one finds that these 
unknown higher order contributions should not hinder the extraction of \V u b\ down to 
an accuracy of 8%. Furthermore we note that the difference between the LO,NLO and 
NNLO curves is minimal for fi m 2.5 GeV, i.e. the apparent convergence of the per- 
turbation series is optimal for this choice of the renormalization scale (which is within 
the above mentioned interval). This fact is not surprising as the size of the coefficients 
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Figure 2: Renormalization scale dependence of T(b — *> X u eV e ) near \x = mb(mb) in 
different orders of perturbative QCD for af\mz) = 0.118, m^\m^) = 4.3 GeV and 
\V ub \ = 0.003. 

hi and h 2 in Eq. ([8]) is greatly reduced for this value of fi. It is also in agreement with 
the observation |3(| that a more natural renormalization scale for the MS renormal- 
ized b quark mass in inclusive semileptonic b decay is somewhat below rrib(mb) as the 
scale is set by the characteristic energy that is released into the hadronic final state. 
Adopting fi »s 2.5 GeV as a choice for the renormalization scale we conclude that the 
a 2 s order approximation to Y{b — > X u eV e ) that is obtained in this letter provides a good 
theoretical foundation for the extraction of \V u b\ from inclusive semileptonic B meson 
decays. 

Useful discussions with K. Chetyrkin, A. Czarnecki, Th. Mannel, K. Melnikov and 
Y.-P. Yao are gratefully acknowledged. This work was supported in part by BMBF 
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